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Abstract. 

A picture for thermodynamics of the glassy state is introduced. It assumes that one 
extra parameter, the effective temperature, is needed to describe the glassy state. This 
explains the classical paradoxes concerning the Ehrenfest relations and the Prigogine- 
Defay ratio. 

As a second part, the approach connects the response of macroscopic observables 
to a field change with their temporal fluctuations, and with the fluctuation-dissipation 
relation, in a generalized non-equilibrium way. 



I INTRODUCTION 

Non-equilibrium thermodynamics for systems far from equilibrium has long been 
a field of confusion. A typical application is window glass. Such a system is far 
from equilibrium: a cubic micron of glass is neither a crystal nor an ordinary under- 
cooled liquid. It is an under-cooled liquid that, in the glass formation process, has 
fallen out of its meta-stable equilibrium. 

Until our recent works on this field, the general consensus reached after more 
than half a century of research was: Thermodynamics does not work for glasses, 
because there is no equilibrium [1] . This conclusion was mainly based on the failure 
to understand the Ehrenfest relations and the related Prigogine-Defay ratio. It 
should be kept in mind that, so far, the approaches leaned very much on equilibrium 
ideas. Well known examples are the 1951 Davies- Jones paper [2], the 1958 Gibbs- 
DiMarzio [3] and the 1965 Adam-Gibbs [4] papers, while a 1981 paper by DiMarzio 
has title "Equilibrium theory of glasses" and a subtitle 'An equilibrium theory of 
glasses is absolutely necessary" [5]. We shall stress that such approaches are not 
applicable, due to the inherent non-equilibrium character of the glassy state. 

Thermodynamics is the most robust field of physics. Its failure to describe the 
glassy state is quite unsatisfactory, since up to 25 decades in time can be involved. 
Naively we expect that each decade has its own dynamics, basically independent 
of the other ones. We have found support for this point in models that can be 
solved exactly. Thermodynamics then means a description of system properties 
under smooth enough non-equilibrium conditions. 



II THERMODYNAMIC PICTURE FOR A SYSTEM 
DESCRIBED BY AN EFFECTIVE TEMPERATURE 



A state that slowly relaxes to equilibrium is characterized by the time elapsed so 
far, sometimes called "age" or "waiting time" . For glassy systems this is of special 
relevance. For experiments on spin glasses it is known that non-trivial cooling or 
heating trajectories can be described by an effective age [6]. Yet we do not wish 
to discuss spin glasses. They have an infinity of long time-scales, or infinite order 
replica symmetry breaking. 

We shall restrict to systems with one diverging time scale, having, in the mean 
field limit, one step of replica symmetry breaking. They are systems with first- 
order-type phase transitions, with discontinuous order parameter, though usually 
there is no latent heat. (As we shall discuss, the same approach applies to true first 
order glassy transitions that do have a latent heat.) 

We shall consider glassy transitions for glass forming liquids as well as for ran- 
dom magnets. The results map onto each other by interchanging volume V, pres- 
sure p, compressibility k = —dlnV/dp, and expansivity a = dlnV/dT, by mag- 
netization M, field H, susceptibility x = {^/N)dM/dH , and "magnetizability" 
a = (-l/N)dM/dT, respectively. 

The picture to be investigated in this work starts by describing a non-equilibrium 
state characterized by three parameters, namely T,p and the age t, or, equivalently, 
T,p and the effective temperature T e (t). This quantity has to follow from solving 
the dynamics of the system, or from doing appropriate experiments. For a set of 
smoothly related cooling experiments Ti(t) at pressures pi, one may express the 
effective temperature as a continuous function: T eji (t) — > T e (T,p). This sets a sur- 
face in (T, T e ,p) space, that becomes multi-valued if one first cools, and then heats. 
For covering the whole space one needs to do many experiments, e.g., at different 
pressures and different cooling rates. The results should agree with findings from 
heating experiments and aging experiments. Thermodynamics amounts to giving 
differential relations between observables at nearby points in this space. 

Of special importance is the thermodynamics of a thermal body at temperature 
T 2 in a heat bath at temperature 7\ = T. This could apply to mundane situations 
such as a cup of coffee, or an ice-cream, in a room. There are also two entropies, 
Si and S 2 - Notice that there are also two time-scales: the time-scale for heat to 
leave the cup is some ten minutes, while the time-scale for equilibrating that heat 
in the room is much smaller. This separation of time-scales allows the difference in 
temperatures. The change in heat of such a system obeys &Q < TxdSx + T 2 dS 2 . 

A similar two-temperature approach proves to be relevant for glassy systems. 
The known exact results on the thermodynamics of systems without currents can 
be summarized by the very same change in heat [7] [8] 



where S ep is the entropy of the fast or equilibrium processes (/3-pro cesses) and I 
the configurational entropy of the slow or configurational processes (a-processes). 



dQ = TdSep + T e dJ 




This object is also known as information entropy or complexity. In the standard 
definition [3] the configurational entropy X is the entropy of the glass minus the one 
of the vibrational modes of the crystal. For polymers it still includes short-distance 
rearrangements, which is a relatively fast mode. It was confirmed numerically 
that X indeed does not vanish at any temperature, thus violating the Adam-Gibbs 
relation between timescale and configurational entropy [9]. In our definition the 
configurational entropy only involves long-time processes; the relatively fast ones 
are counted in S ep . With this point of view the applicability of the Adam-Gibbs 
relation remains an open issue. 

It is both surprising and satisfactory that a glass can be described by the same 
general law. If, in certain systems, also an effective pressure or field would be 
needed, then 6.Q is expected to keep the same form, but dW would change from 
its standard value —pdV for liquids, or —MdH for magnets. In the latter case 
it would become — Mi&H — M 2 dH e , where H e is the effective field, and M 1 and 
M 2 add up to M. Such an extension could be needed to describe a larger class of 
systems. 



A First and second law 

For a glass forming liquid the first law dU = dQ + dW becomes 

dU = TdS cp + T e dX - pdV 
It is appropriate to define the generalized free enthalpy 



G = U - TS, 



cp 



T e T + pV 



This is not the standard form, since T e ^ T. It satisfies 

dG = -S cp dT - XdT e + Vdp 

The total entropy is 

S = S^p + X 

The second law requires dQ < TdS, leading to 

(T e - T)dX < 0, 



(2) 



(3) 



(4) 



(5) 



(6) 



which merely says that heat goes from high to low temperatures. 

Since T e = T e (T,p), and both entropies are functions of T, T e and p, the expres- 
sion (1) yields a specific heat 
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In the glass transition region all factors, except dxT e) are basically constant. This 
leads to 



C P = C 1 + C 2 ^\ p (8) 

Precisely this form has been assumed half a century ago by Tool [10] as starting 
point for the study of caloric behavior in the glass formation region, and has often 
been used for the explanation of experiments [2] [11]. It is a direct consequence of 
eq. (1). 

For magnetic systems the first law brings 

dU = TdS cp + T e dX - MdH (9) 

As above, one can define the free energy F = U — TS cp — T e X. It satisfies the 
relation dF = -S cp dT - XdT e - MdH. 



B Modified Maxwell relation 

For a smooth sequence of cooling procedures of a glassy liquid, eq. (2) implies 
a modified Maxwell relation between macroscopic observables such as U(t,p) — > 
U(T,p) = U(T,T e (T,p),p) and V. This solely occurs since T e is a non-trivial 
function of T and p for the smooth set of experiments under consideration. 

For glass forming liquids it reads 
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This is the modified Maxwell relation between observables U and V. In equilibrium 
T e = T, so the right hand side vanishes, and the standard form is recovered. 
Similarly, one finds for a glassy magnet 
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C Modified Clausius-Clapeyron relation 

Let us consider a first order transition between two glassy phases A and B. An 
example could be the transition from low-density-amorphous ice to high-density- 
amorphous ice [12]. For the standard Clausius-Clapeyron relation one uses that the 
free enthalpy G is continuous along the first order phase transition line p g {T). Since 
T e 7^ T, it is actually not obvious that the function G should still be continuous 
there. Based on a quasi-static approach involving a partition sum and, equivalently, 
on the experience that in mean field models replica theory has always brought the 



relevant physical free energy, we expect that our generalized free enthalpy (3) is 
indeed continuous at this transition. 

Let us consider a first order transition between phases A and B. Assuming that 
AG = along the glass transition line p g (T), a little amount of algebra yields the 
modified Clausius-Clapeyron relation 

AT dp q AU+p q AV . (dT e ^ T e \ 

dT T \dT T ) y ' 

where d/dT = d/dT + (dp g /dT)d/dp is the "total" derivative along the transition 
line. 

It would be very interesting to test this relation for the two forms of amorphous 
ice. For those substances Mishima and Stanley [13] have presented a thermody- 
namic construction of the standard free enthalpy or Gibbs potential G. It is, how- 
ever, based on equilibrium ideas and, in contrast to eq. (3), does not involve the 
effective temperature in the amorphous phases. This approach therefore predicts 
the validity of the standard Clausius-Clapeyron relation. We feel that the standard 
G is not the physically relevant one, and that the analysis should be redone, by 
considering many cooling rates and taking into account the effective temperature 
and the possible violation of the Clausius-Clapeyron relation by going from the 
high density amorphous phase to the low density one, and vice versa. 

When phase A is an equilibrium under-cooled liquid, and phase B is a glass, it 
holds that T e = T in phase A. Then the relation reduces to 

{Va-Vb)-^ = j, h (— - — Ji-B (13) 

Several models for magnets undergoing a first order glassy transitions with a latent 
heat have been studied in the literature [14] [15] [16] [17]. 



D Ehrenfest relations and Prigogine-Defay ratio 

In the glass transition region a glass forming liquid exhibits smeared jumps in the 
specific heat C p , the expansivity a and the compressibility k. If one forgets about 
the smearing, one may consider them as true discontinuities, yielding an analogy 
with continuous phase transitions of the classical type. 

Following Ehrenfest one may take the derivative of AV(T, p g (T)) = 0. The result 
for a glass forming liquid may be written as 

Aa = Ak^ (14) 

while for a glassy magnet 



(15) 



The conclusion drawn from half a century of research on glass forming liquids is 
that this relation is never satisfied [2] [18] [19] [1]. This has very much hindered 
progress on a thermodynamical approach. However, from a theoretical viewpoint 
it is hard to imagine that something could go wrong when just taking a derivative. 
We have pointed out that this relation is indeed satisfied automatically [7], but it 
is important say what is meant by k in the glassy state. 

Let us make an analogy with spin glasses. In mean field theory they have infinite 
order replica symmetry breaking. From the early measurements of Canella and 
Mydosh [20] on AuFe it is known that the susceptibility depends logarithmically 
on the frequency, so on the time scale. The short-time value, called Zero-Field- 
Cooled (ZFC) susceptibility is a lower bound, while the long time value, called 
Field-Cooled (FC) susceptibility is an upper bound. Let us use the term "glassy 
magnets" for spin glasses with one step of replica symmetry breaking. They are 
relevant for comparison with glass forming liquids. For them the situation is worse, 
as the ZFC value is discontinuous immediately below T g . This explains why already 
directly below the glass transition different measurements yield different values for 
k. These notions are displayed in figure 1. 
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FIGURE 1. Schematic plot of the field-cooled (FC) and zero-field-cooled (ZFC) susceptibility 
in realistic spin glasses and in glassy magnets, as function of temperature, in arbitrary units. 
In realistic spin glasses the infinite time or field-cooled susceptibility is larger than the short 
time or zero-field-cooled susceptibility. In magnetic analogs of realistic glasses the short time 
susceptibility even has a smeared discontinuity at the glass transition, yielding a value of \ that 
depends on the precise type of experiment which is performed. In glass forming liquids the same 
happens for the compressibility. 



Previous claims about the violation of the first Ehrenfest relation can be traced 
back to the equilibrium thermodynamics idea that there is one, ideal k, to be 
inserted in (14). Indeed, investigators always considered cooling curves V(T,pi) 
at a set of pressures pi to determine Aa and dp g /dT. However, An was always 
determined in another way, often from measurements of the speed of sound, or 
by making more complicated pressure steps [21]. In equilibrium such alternative 
determinations would yield the same outcome. In glasses this is not the case: 
the speed of sound is a short-time process, and additional pressure steps modify 
the glassy state. Therefore alternative procedures should be avoided, and only the 
cooling curves V(T,pi) should be used. They constitute a liquid surface V\ iqui d(T, p) 
and a glass surface V g i ass (T, p) in (T,p,V) space. These surfaces intersect, and 
the first Ehrenfest relation is no more than a mathematical identity about the 
intersection line of these surfaces. It is therefore automatically satisfied [7]. The 
most careful data we came across were collected by Rehage and Oels for atactic 
polystyrene [21]. In figure 2 we present those data in a 3-d plot, underlining our 
point of view. 

The second Ehrenfest relation follows from differentiating AU(T,p g (T)) = 0. 
The obtained relation will also be satisfied automatically. However, one then elim- 
inates dU /dp by means of the Maxwell relation. We have already discussed that 
outside equilibrium it is modified. As a result we obtain 

AC P dp 1 / 8T e \ dl 

T g V =Aa dT + v[ 1 -W P ) dT (16) 

The X term constitutes the "total" derivative of the configurational entropy along 
the glass transition line. Its prefactor only vanishes at equilibrium, in which case 
the standard Ehrenfest relation is recovered. 
For glassy magnets one has similarly 
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The equality T e (T,p g (T)) = T implies the resulting identity 



(17) 



dT e dT e 



dT dT 



dT 



dp 



d Pg 



T dT y ' 



Combining the two Ehrenfest relations one may eliminate the slope of the transition 
line. This leads to the so-called Prigogine-Defay ratio 

n= ACpAK (19) 

TV(Aa) 2 [ } 



This looks like an equilibrium quantity. For equilibrium transitions it should be 
equal to unity. Assuming that at the glass transition a number of unspecified 



parameters undergo a phase transition, Davies and Jones derived that II > 1 [2], 
while DiMarzio showed that in that case the correct value is II = 1 [22]. In glasses 
typical experimental values are reported in the range 2 < II < 5. It was therefore 
generally expected that II > 1 is a strict inequality arising from the requirement of 
mechanical stability. 

We have pointed out, however, that, since the first Ehrenfest relation is satisfied, 
it holds that 



AC P 1 / dT e 
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(20) 



Depending on the smooth set of experiments to be performed, dp g /dT can be small 
or large: II depends on the set of experiments. As a result, it can also be below 
unity. Rehage-Oels found II = 1.09 ~ 1 at p — 1 kbar, using a short-time value for 
k. Reanalyzing their data we find from (20), where the correct n has been inserted, 




FIGURE 2. Data of the glass transition for cooling atactic polystyrene at rate 18 K/h, scanned 
from the paper of Rchage and Oels (1976): specific volume V (cm 3 /g) versus temperature T (K) 
at various pressures p (kbar). As confirmed by a polynomial fit, the data in the liquid essentially 
lie on a smooth surface, and so do the data in the glass. The first Ehrenfest relation describes no 
more than the intersection of these surfaces, and is therefore automatically satisfied. The values 
for the compressibility derived in this manner will generally differ from results obtained via other 
procedures. 



a value II = 0.77, which indeed is below unity. The commonly accepted inequality 
II > 1 is based on the equilibrium assumption of a unique k. Our theoretical 
arguments and the Rehage-Oels data show that such idea's are incorrect. 



E Fluctuation formula 



The basic result of statistical physics is that it relates fluctuations in macroscopic 
variables to response of their averages to changes in external field or temperature. 
We have wondered whether such relations generalize to the glassy state. We have 
found arguments in favor of such a possibility both from the fluctuation-dissipation 
relation and by exactly solving the dynamics of model systems [23]. Susceptibilities 
appear to have a non-trivial decomposition, that looks as being very general. Here 
we give arguments leading to it. 

In cooling experiments at fixed field it holds that M =M(T(t),T e (t, H), H). For 
thermodynamics one eliminates time, implying M =M(T,T e (T, H), H). One may 
then expect two terms: 
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The first is the fluctuation contribution 
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To calculate it, we switch from a cooling experiment to an aging experiment at 
the considered T, T e and H, by keeping, in Gedanken, T fixed from then on. The 
system will continue to age, expressed by T e = T e (t;T,H). We may then use the 
equality 
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We have conjectured [23] that the left hand side may be written as the sum of 
fluctuation terms for fast and slow processes, 
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(24) 



The first term is just the standard equilibrium expression for the fast equilibrium 
processes. Notice that the slow processes enter with their own temperature, the 
effective temperature. This decomposition is confirmed by use of the fluctuation- 
dissipation relation in the form to be discussed below. Combination of (22), (23) 
and (24) yields our non-equilibrium prediction 
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The first two terms are instantaneous, and thus the same for aging and cooling. 
The third term is a correction, related to an aging experiment. In the models 
considered so far, this term is small [23] [24]. 

Since T e 7^ T, there occurs in eq. (21) also a new, configurational term 



x conf = 
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It originates from the difference in the system's structure for cooling experiments at 
nearby fields. This is the term that is responsible for the discontinuity of \ & t the 
glass transition. For glass forming liquids such a term occurs in the compressibility. 
Its existence was anticipated in some earlier work [18] [19]. 



F Fluctuation-dissipation relation 

Nowadays quite some attention is payed to the fluctuation-dissipation relation in 
the aging regime of glassy systems. It was first put forward in works by Horner [25] 
[26] and generalized by Cugliandolo and Kurchan, see [27] for a review. 

In the aging regime there holds a fluctuation-dissipation relation between the cor- 
relation function C(t,t')=< 5M(t)5M(t') > and G(t,t'), the response of < M(t) > 
to a short, small field change 5H(t') applied at an earlier time 

—^r^~ = T e (t,t')G(t,t') (27) 

with T e (t,t') being an effective temperature, formally defined by this relation. In 
the equilibrium or short-time regime \t — t'\ <C t, it is just equal to T; in the aging 
regime t/t' = 0(1), it depends only logarithmically on t and t', making it a useful 
concept. 

We have observed that in simple models without fast processes T e (t,t f ) = T e (t') 
is a function of one of the times only. One then expects that T e (t) is close to the 
"thermodynamic" effective temperature T e (t). We have shown that [24] 

T e (t) = T e (t)-T e (t)i^ g^l =t ) +■■■ (28) 

So the effective temperatures T e and T e are not identical. However, in the models 
analyzed so far, the difference is subleading in 1/lnt. 

Notice that the ratio d t iC(t,t')/G(t,t f ) = T e (t') is allowed to depend on time 
t'. The situation with constant T e is well known from mean field spin glasses [27], 
but we have not found such a constant T e beyond mean- field [23] [24]. Only at 
exponential time-scales the mean field spin glass behaves as a realistic system [8]. 



G Time-scale arguments 



Consider a simple system that has only one type of processes (a processes), which 
falls out of equilibrium at some low T. When it ages a time t at T — it will have 
achieved a state with effective temperature T e , that can be estimated by equating 
time with the equilibrium time-scale. Let us define T e by 

t = r eq {T e ) (29) 

We have checked in solvable models that, to leading order in mi, it holds that 
T e = T e . (The first non-leading order turns out to be non- universal). This equality 
also is found in cooling trajectories, when the system is well inside the glassy 
regime. It says that the system basically has forgotten its history, and ages on its 
own, without caring about the actual temperature. Another way of saying is that 
dynamics in each new time-decade is basically independent of previous decade. 

In less trivial systems, for instance those having a Vogel-Tammann-Fulcher law, 
the time-scale may have parameters that depend on the actual temperature, im- 
plying r = r(T, T e ). We have already found support for the expectation that, to 
leading order, T e follows by equating this expression with time t. 

In many systems one finds a t'/t scaling in the aging regime of two-time quantities. 
There is a handwaving argument to explain that: 

c <*^ MC <;^> w0 <^> = <?> (30) 

showing indeed the familiar t/t' scaling. In the models studied so far we have found 
logarithmic scaling corrections [23] [24]. 

Ill SOLVABLE MODELS 

The above relations have been tested in models of which the statics is exactly 
solvable, and the dynamics is partially solvable, such as the p-spin model [26] [7] 
and a directed polymer model with glassy behavior [28] [7]. Partial results follow 
from the backgammon model [29], [30]. 

More interesting are models with exactly solvable parallel Monte Carlo dynamics, 
such as independent harmonic oscillators [31] [23] [24] or independent spherical spins 
in a random field [23] [24]. In both cases a glassy behavior occurs when cooling 
towards low temperatures. A related model with a set of fast modes and a set of 
slow modes, that have a Vogel-Fulcher-Tammann-Hesse law for the divergence of 
the equilibrium time-scale, is currently under study. 

The two-temperature approach put forward here also explains the thermodynam- 
ics of black holes [32] and star clusters [33]. 
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